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Abstract 

We calculate electron-hole generation and recombination rates for Coulomb scattering (Auger 
recombination and impact ionization) in Graphene. The conduction and valence band dispersion 
relation in Graphene together with energy and momentum conservation requirements restrict 
the phase space for Coulomb scattering so that electron-hole recombination times can be much 
longer than 1 ps for electron- hole densities smaller than 10 12 cnr 2 . 
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1 Introduction 



Graphene is a single two dimensional (2D) atomic layer of carbon atoms forming a dense 
honeycomb crystal lattice [1J . The electronic properties of Graphene have generated tremendous 
interest in both experimental and theoretical arenas [21 El El El E]- The massless energy 
dispersion relation of electrons and holes with zero (or close to zero) bandgap results in novel 
behavior of both single-particle and collective excitations [H [21 El El El El E] • The high mobility of 
electrons in Graphene has prompted theoretical and experimental investigations into Graphene 
based ultra high speed electronic devices such as field-effect transistors, pn-junction diodes, and 
terahertz oscillators O El El [El CEH IE] • The behavior of many of these devices depends on the 
electron-hole recombination rates in Graphene. For example, the diffusion length of injected 
minority carriers in a pn-junction diode is proportional to the square-root of the minority 
carrier recombination time [13] . It is therefore important to understand the mechanisms that 
are responsible for electron-hole generation and recombination in Graphene and the associated 
time scales. 

Small band-gap semiconductors usually have large electron-hole recombination rates due 
to Coulomb scattering (Auger recombination) [TJ]. Graphene, with a zero (or close to zero) 
bandgap, presents a limiting case. The zero bandgap and the large optical phonon energy in 
Graphene (196 meV [15] ) suggest that electron-hole recombination rates could be dominated 
by Auger processes. In addition, the zero bandgap also implies that electron-hole generation 
rates in Graphene due to Coulomb scattering (impact ionization) may also not be small even 
in the absence of high energy carriers. 

In this paper we calculate the electron-hole generation and recombination rates for Coulomb 
scattering (Auger recombination and impact ionization) in Graphene. We show that the 
conduction and valence band dispersion relation in Graphene together with energy and 
momentum conservation requirements restrict the phase space for Coulomb scattering so that 
electron-hole generation recombination times can be much longer than 1 ps at all temperatures 
for electron-hole densities smaller than 10 12 cm -2 . 
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band 

Figure 1: Electron-hole recombination in Graphene from Coulomb scattering (Auger 
recombination) via the CCCV and the CVVV processes is shown. The two processes shown 
are mirror images of each other. 




band 

Figure 2: Electron- hole generation in Graphene from Coulomb scattering (impact ionization) 
via the CVCC and the VVCV processes is shown. The two processes shown are mirror images 
of each other. In addition, the generation processes shown are the reverse of the recombination 
processes shown in Fig.l. 
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2 Theoretical Model 



In Graphene, the valence and conduction bands resulting from the mixing of the p 2 -orbitals are 
degenerate at the inequivalent K and K' points of the Brillouin zone pQ. Near these points, 
the conduction and valence band dispersion relations can be written compactly as, 

E s (k) = shv\k\ (1) 

where s = ±1 stand for conduction (+1) and valence (—1) bands, respectively, and v is the 
"light" velocity of the massless electrons and holes, and assumed to be equal to 10 8 cm/s. The 
wavevector k is measured from the K(K') point. 

Electron-hole recombination due to Coulomb scattering (Auger recombination) in Graphene 
occurs by the two processes depicted in Fig.l. In the CCCV process (Fig.l), an electron in the 
conduction band with initial momentum ki scatters off another electron in the conduction band 
with momentum k2. The result is an electron in the conduction band with momentum ki + Q 
and an electron in the valence band with momentum k2 — Q. In the CVVV process (Fig.l), a 
hole in the valence band with initial momentum ki scatters off another hole in the valence band 
with momentum k2 . The result is a hole in the valence band with momentum ki + Q and a hole 
in the conduction band with momentum k2 — Q. The CVVV process is a mirror image of the 
CCCV process. The labeling scheme used here indicates the initial and final states of the two 
electrons involved in the scattering process. The electron- hole recombination rates RcccY( ri ',p) 
and Rcwv(n,p) are functions of the electron and hole densities, n and p, respectively, and 
the symmetry between the conduction and valence band dispersions in Graphene implies that 
Rcccvi^iP) = Rcvvv(p,n). Electron- hole generation due to Coulomb scattering (impact 
ionization) occurs by the reverse of the CCCV and CVVV recombination processes. We label 
these generation processes as CVCC and VVCV, respectively, and they are depicted in Fig. 2. 
The generation rates, Gcvcc( n ,p) an d GVvcv(p, n ), also satisfy Gcvcc( n ,p) = Gyycv(p,n). 
In thermal equilibrium, the total generation and recombination rates must be equal, i.e. 
Gcvcc + GWcv = -Rcccv + -Rcvvv- 
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Given the symmetry between the various processes, it is enough to consider in detail only 
one process. In the discussion that follows we will concentrate only on the CCCV Auger 
recombination process. In the CCCV process, energy conservation implies, 

fru|ki| + fru|k 2 | = hv\ki + Q| - Hv\k 2 - Q| 

|ki| + |k 2 | = |ki + Q| - |k 2 -Q| (2) 

The available phase space for Auger recombination can be understood as follows. For any three 
vectors ki, k 2 , and Q one has the identity, 

|ki + Q| - |k 2 - Q| < |ki + k 2 | < |ki| + |k 2 | (3) 

The energy conservation condition requires both the inequalities above to be equalities. The 
inequality on the right would be an equality if (and only if) the vectors ki and k 2 point in 
the same direction. The inequality on the left would be an equality if (and only if) the vectors 
ki + Q and k 2 — Q point in the opposite direction. If ki and k 2 point in the same direction, 
then ki + Q and k 2 — Q will point in the opposite direction if (and only if) ki and Q also 
point in the same direction and |Q| > |k 2 |. Energy conservation therefore requires that the 
vectors ki, k 2 , and Q all lie on the same line. This requirement also holds for all the Coulomb 
scattering processes depicted in Fig.l and Fig. 2. 

The Bloch functions for the conduction (s = +1) and valence (s = —1) band electrons in 
Graphene can be written as pQ, 

gik.r 

tp s ,k( r ) = u S}k (r) (4) 

Here, N is the total number of unit cells in the crystal. The periodic part u s ^(r) of the Bloch 
function has the following overlap integral [I], 

I < iV,k'K,k > | 2 = - [1 + as' cos(0 k ',k)] (5) 

where #k',k is the angle between the vectors k' and k. We assume that the occupation statistics 
of electrons in conduction and valence bands are given by the Fermi-Dirac distribution functions 
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/ s (k), 



/s(k) ~ 1 + e (E a (k)-E fs )/KT ( 6 ) 



Ef s are the Fermi levels. We assume different Fermi levels for conduction and valence electrons 
to allow for non-equilibrium electron-hole populations, as is the case in a forward biased 
pn-junction diode [13J. The electron and hole densities are given as follows, 



n = 4 /(^ /+l(k) (7) 



P 



'/csp 11 -'- 10 " 1 (8) 



The factor of 4 in the front accounts for spin degeneracy and the two valleys at K and K' . 

The electron-hole recombination rate Rcccv{"n,p) (units: cm _2 -sec _1 ) due to Auger 
scattering can be written as [17J, 

WO,,) = 2 (^)/|^/0/0,M( kl , k2 , Q ,p 

[1 - /_!(k 2 - Q)] [1 - / +1 (k! + Q)] / +1 (k0 /+i(k 2 ) 

8(hv\k. x \ + hv\k 2 \ - hv\k! + Q| + hv\k 2 - Q|) (9) 

The factor of two in the front comes from the two valleys at K and K' . The scattering matrix 
element M(ki,k 2 , Q) includes both direct and exchange processes, and can be written as, 

|M(ki, k 2 , Q)| 2 = |Afd(ki,k 2 , Q)| 2 + |M e (ki, k 2 , Q)| 2 + |M d (k x ,k 2 , Q) - M e (ki, k 2 , Q)| 2 (10) 

Assuming statically screened Coulomb interaction, the matrix elements, M<j(ki,k 2 ,Q) and 
M e (ki,k 2 ,Q), for the direct and exchange scattering processes, respectively, are as follows, 

e 2 

M d (ki,k 2 ,Q) = — - < u+i,ki+q|«+i,ki >< u-i,k2-Q|w+i,k 2 > ( n ) 

^£oo {\m + qtf) 

e 2 

M e (k 1)k2 ,Q) = 2£oo (|Q + kl -k 2 |+Q TF ) < u +lM+Q\u + lte >< U- 1M - Q \u +1M > 

(12) 

Here, e is the electron charge, Qtf is the Thomas-Fermi wavevector |18j . and eoo is the average 
of the dielectric constants of the media on both sides of the Graphene layer. The relative 
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directions for the vectors ki, k2, and Q allowed by energy conservation results in the values of 
all the overlap integrals in Equations (jlip and (|12p to equal unity. Assuming screening by both 
electrons and holes, the expression for the Thomas- Fermi wavevector in Graphene is [18], 

Qtf = ^r~ 2 log \(e E ^' KT + l) L- E ^' KT + 1)1 (13) 

where, -Ey+i and are the Fermi levels for the conduction and valence electrons, respectively. 
After integrating out the delta function, the six-dimensional integral in Equation can be 
reduced to the following three-dimensional integral, 

1 roo Ah, roo Ah roo Af) , 

Rcccv(n,p) = T2- ~^T~ / ^ \M(h,k 2 , Q)| 2 J(h + Q)(Q- k 2 ) h k 2 

fi v Jo ^tt Jo to Jk 2 to 

[i - f-i(Q - h)) [1 - f +1 (h + Q)} f +1 (h) f +1 (k 2 ) (14) 

The above equation is the main result of this paper. The total Auger recombination rate R(n,p) 
is the sum of the rates of the CCCV and CVVV processes, 

R(n,p) = Rcccv(n,p) + R C vvv(n,p) (15) 



The average electron-hole recombination time r r due to Coulomb scattering is defined as, 

1 R(n,p) 



(16) 



T r min(n,p) 

where the smaller carrier density appears in the denominator on the right hand side. r r can 
also be interpreted as the minority carrier lifetime in situations where the electron and hole 
densities are very different (see the discussion below). Using the result in Equation (|14p . the 
generation rate Gcvcc(n,p) due to CVCC process can be written as, 

1 roo Ah, roo Ah roo Af) , 

G C vcc(n,p) = -s- / -± \M(h,k 2 , Q)\ 2 J(h + Q) (Q - k 2 ) k x k 2 

Tl V JO to Jo Z7T Jk 2 to v 

[1 - f + i(h)] [l - f +1 (k 2) ] /_i(Q - k 2 ) f +1 (Q + kt) (17) 
The total generation rate G(n,p) is the sum of the rates of the CVCC and VVCV processes, 

G(n,p) = G C vcc(n,p) + Gwcv(n,p) (18) 
The average electron-hole generation time r g due to Coulomb scattering is defined as, 

I = °M) (19) 

T g mm{n,p) 
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Figure 3: Electron-hole generation (dashed) and recombination (solid) rates are plotted as a 
function of the electron and hole densities (assumed to be equal) for different temperatures. 
The different curves correspond to temperatures T=10K, 77K, and 300K. The assumed values 
of v and are 10 s cm/s and 10e o , respectively. 




Figure 4: Electron-hole generation (dashed) and recombination (solid) times are plotted as a 
function of the electron and hole densities (assumed to be equal) for different temperatures. 
The different curves correspond to temperatures T=10K, 77K, and 300K. The assumed values 
of v and are 10 8 cm/s and 10e o , respectively. 
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Figure 5: Electron-hole generation (dashed) and recombination (solid) times are plotted as a 
function of the electron and hole densities (assumed to be equal) for different temperatures. 
The different curves correspond to temperatures T=10K, 77K, and 300K. The assumed values 
of v and Eqo are 10 s cm/s and 4e G , respectively. 

3 Results 

Fig. 3 shows the total generation and recombination rates plotted as a function of the electron 
and hole densities, which are assumed to be equal, for different temperatures (T=10K, 77K, 
and 300K). The value of eoo used in simulations is 10e o assuming Aluminum-oxide on both 
sides of the Graphene layer [9]. Fig. 4 shows the corresponding generation and recombination 
times. For any temperature, the generation and recombination rate curves cross at the point 
where the electron and hole densities have their thermal equilibrium values. Just like in two 
dimensional semiconductor quantum wells, the temperature dependence of Coulomb scattering 
is a sensitive function of the electron and hole densities as well as temperature |14| . At higher 
temperatures the probability of finding energetic electrons and holes is larger and therefore 
the generation rate increases with temperature. For small electron-hole densities, increase in 
temperature spreads the carrier distributions to higher energies where Auger recombination 
is less efficient and therefore the recombination rate decreases. But for large electron-hole 
densities, the electrons and holes near the band edges can recombine only if the final scattering 
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Figure 6: Electron-hole generation (dashed) and recombination (solid) rates are plotted for 
different electron densities. The hole density is fixed and equals 10 12 cm -2 . T=300K. The 
contributions Gcvcc and GVvcv to the total generation rate G, and the contributions i?cccv 
and i?cvw to the total recombination rate R, are also plotted. 

states are unoccupied. An increase in temperature generates more unoccupied states. As a 
result of the above two factors, the recombination rates at large electron-hole densities are less 
sensitive to temperature. 

Fig. 4 shows that for electron-hole densities smaller than 10 12 cm -2 the recombination time 
is longer than 1 ps at all temperatures, and for electron-hole densities smaller than 10 11 cm -2 
the recombination time is longer than 5 ps at all temperatures. Fig. 5 shows the generation 
and recombination times for = 4e Q assuming Silicon-dioxide on both sides of the Graphene 
layer. Clearly, the Coulomb scattering rates depend on the dielectric surrounding the Graphene 
layer and are enhanced for a smaller dielectric constant medium. The effect of the surrounding 
dielectric on the generation and recombination times is more pronounced at small electron-hole 
densities when carrier screening of the Coulomb interaction is less effective. 

Fig.6 shows the individual contributions, i?cccv and -Rcvvv of the CCCV and CVVV 
processes, respectively, to the total Auger recombination rate R, and the contributions Gcvcc 
and Gwcv of the CVCC and VVCV processes, respectively, to the total generation rate G. The 
electron density is varied and the hole density is fixed at 10 12 cm~ 2 . T=300K and eoo = 10e o . 
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Figure 7: Minority carrier (electron) lifetimes is plotted as a function of the majority (hole) 
density for different temperatures (T=10K, 77K, and 300K). The assumed values of v and eoo 
are 10 8 cm/s and 10e o , respectively. 

The thermal equilibrium value of the electron density, corresponding to the hole density of 10 12 
cm~ 2 , is 1.5 x 10 9 cm -2 . For electron densities much smaller than the hole density the CVVV 
process dominates recombination, and the recombination rate varies approximately linearly 
with the electron density. For electron densities much larger than the hole density the CCCV 
process dominates recombination. When electron and hole densities are equal then, as explained 
earlier, Rqccv = ^CWV- For electron densities much smaller than the hole density the VVCV 
process dominates generation and the generation rate is almost independent of the electron 
density. For electron densities much larger than the hole density the CVCC process dominates 
generation. For equal electron and hole densities, Gcvcc = GVvcv- 

For device applications, it is also interesting to look at the minority carrier 
generation-recombination rates in situations where the electron (or the hole) density is much 
smaller than the hole (or the electron) density. This is the case, for example, in a forward biased 
pn-junction [13]. Fig. 6 and the discussion above shows that the minority carrier generation and 
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recombination rates can be written approximately as follows, 

R - G = ^—^ (20) 

where, we have assumed that the electrons are the minority carriers and holes are the majority 
carriers. n Q is the thermal equilibrium electron density and r r is the minority carrier (electron) 
lifetime. The minority carrier lifetime is independent of the minority carrier density but depends 
on the majority carrier density and the temperature. From Fig. 6, the minority carrier lifetime 
for a majority carrier density of 10 12 cm -2 is approximately 1.1 ps at T=300K. Fig.7 shows 
the minority carrier (electron) lifetime as a function of the majority carrier (hole) density for a 
different temperatures (T=10K, 77K, and 300K). eoo = 10e o . 

In conclusion, we have calculated electron-hole generation and recombination rates due to 
Coulomb scattering in Graphene. Our results show that electron-hole recombination times in 
Graphene can be much longer than 1 ps at all temperatures for electron-hole densities smaller 
than 10 12 cm" 2 . The author would like to acknowledge helpful discussions with Edwin Kan 
and Sandip Tiwari. 
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